In this notebook we will add TP53 alteration status as discussed in #837:
TP53 altered - loss, if :
- A sample contains a TP53 hotspot SNV mutation. (Cancer hotspot database and downloadable file available). Please also crosscheck that all mutations from this table are included.
- A sample contains two TP53 alterations, suggesting (but not confirming) that both alleles are affected (SNV+SNV, (CNV or SV) + SNV).
- A sample contains one alteration (SNV or (CNV or SV)) + has cancer_predispositions == “Li-Fraumeni syndrome”, suggesting there is a germline variant in addition to the somatic variant we observe.
- A sample does not have a TP53 alterations, but has cancer_predispositions == “Li-Fraumeni syndrome” and TP53 classifier score for matched RNA-Seq > 0.5 (or higher cutoff we decide upon later).
Note: CNV and SV will be considered as the same event, but we will be adding SV_counts and SV.type to the altered status output file
TP53 altered - activated, if:
- A sample contains one of the two TP53 activating mutations p.R273C and p.R248W. [@doi:10.1038/ng0593-42]
Setup
library("ggpubr")
Loading required package: ggplot2
Loading required package: magrittr
library("ggthemes")
library("tidyverse")
── Attaching packages ────────────────────────────────── tidyverse 1.2.1 ──
✔ tibble 2.1.3 ✔ purrr 0.3.2
✔ tidyr 0.8.3 ✔ dplyr 0.8.3
✔ readr 1.3.1 ✔ stringr 1.4.0
✔ tibble 2.1.3 ✔ forcats 0.4.0
── Conflicts ───────────────────────────────────── tidyverse_conflicts() ──
✖ tidyr::extract() masks magrittr::extract()
✖ dplyr::filter() masks stats::filter()
✖ dplyr::lag() masks stats::lag()
✖ purrr::set_names() masks magrittr::set_names()
library("ggfortify")
library("broom")
# rootdir
root_dir <- rprojroot::find_root(rprojroot::has_dir(".git"))
data_dir <- file.path(root_dir, "data")
input_dir <- file.path(root_dir,
"analyses",
"tp53_nf1_score",
"input")
# cell composition
cell_line_composition <- read_tsv(file.path("..",
"molecular-subtyping-HGG",
"input",
"cell-line-composition.tsv"),
col_types =
readr::cols( aliquot_id = readr::col_character()))
# histology
if ( params$base_run ==0 ){
clinical<-read.delim(file.path(data_dir,"histologies.tsv"), stringsAsFactors = FALSE)
} else{
clinical<-read.delim(file.path(data_dir,"histologies-base.tsv"), stringsAsFactors = FALSE)
}
histology <- clinical %>%
select("Kids_First_Biospecimen_ID",
"sample_id",
"cancer_predispositions",
"sample_type",
"experimental_strategy",
"composition") %>%
# merge cell line composition
left_join(cell_line_composition)
Joining, by = c("Kids_First_Biospecimen_ID", "sample_id")
# Cancer database
hotspot_database_2017_tp53_snv <- readxl::read_xls(file.path(input_dir,"hotspots_v2.xls"),sheet = 1) %>%
filter(Hugo_Symbol == "TP53")
hotspot_database_2017_tp53_indel <- readxl::read_xls(file.path(input_dir,"hotspots_v2.xls"),sheet = 2) %>%
filter(Hugo_Symbol == "TP53")
# p53, p63 and p73 functional sites
functional_sites_tp53 <- read_tsv(file.path(input_dir,"hotspot_Chen_2006.tsv"))
Parsed with column specification:
cols(
p53 = col_character(),
p63 = col_character(),
p73 = col_character()
)
results_dir <- file.path(root_dir,
"analyses",
"tp53_nf1_score",
"results")
if (!dir.exists(results_dir)) {
dir.create(results_dir)
}
# classifier scores
classifier_score <- read_tsv(
file.path(
results_dir,
"gene-expression-rsem-tpm-collapsed_classifier_scores.tsv"))%>%
select(sample_id,tp53_score) %>%
dplyr::rename(Kids_First_Biospecimen_ID_RNA=sample_id) %>%
left_join(histology ,by=c("Kids_First_Biospecimen_ID_RNA"="Kids_First_Biospecimen_ID"))
Parsed with column specification:
cols(
sample_id = col_character(),
ras_score = col_double(),
tp53_score = col_double(),
nf1_score = col_double(),
ras_shuffle = col_double(),
tp53_shuffle = col_double(),
nf1_shuffle = col_double()
)
# Copy number per sample
# overlapping functional domain
cnv_domain_overlap <- read_tsv(
file.path(
results_dir,
"loss_overlap_domains_tp53.tsv"))
Parsed with column specification:
cols(
biospecimen_id = col_character(),
copy_number = col_double(),
ploidy = col_double(),
domain = col_character(),
sample_id = col_character()
)
# Structural variant overlapping
# or within gene locus of TP53
sv_overlap <- read_tsv(
file.path(
results_dir,
"sv_overlap_tp53.tsv"))
Parsed with column specification:
cols(
SV.chrom = col_double(),
SV.start = col_double(),
SV.end = col_double(),
SV.length = col_double(),
Kids.First.Biospecimen.ID.Tumor = col_character(),
SV.type = col_character(),
Gene.name = col_character(),
ALT = col_character(),
sample_id = col_character()
)
Check if all function sites are in hotspots
All functional sites are just 1 base so checking in hotspot_database_2017_snv
functional_sites_tp53 %>%
filter(!gsub("[A-Z|a-z]","",p53) %in% hotspot_database_2017_tp53_snv$Amino_Acid_Position )
2 functional sites are missing in hotspots database.
SNV for TP53
Removing Silent or Intron classified variants to capture only putative damaging mutations
consensus_tp53_snv_indel <- data.table::fread(
file.path(data_dir,"snv-consensus-plus-hotspots.maf.tsv.gz"),
select = c("Chromosome",
"Start_Position",
"End_Position",
"Strand",
"Variant_Classification",
"Tumor_Sample_Barcode",
"Hugo_Symbol",
"HGVSp_Short"),
data.table = FALSE) %>%
filter(Hugo_Symbol == "TP53") %>%
filter(!(Variant_Classification %in% c("Silent", "Intron",
# remove other non-amino acid SNVs
"3'Flank" ,"5'Flank",
"3'UTR", "5'UTR" )))
Registered S3 method overwritten by 'R.oo':
method from
throw.default R.methodsS3
Gather annotation for SNV
hotspots : overlaps AA position which are statistically significant SNV/Indels activating : overlaps AA position which are found to act as gain-of-function according to literature
consensus_tp53_snv_indel <- consensus_tp53_snv_indel %>%
mutate(
hotspot = case_when(
# strip REF and Variant AA to get Amino_Acid_Position in consensus maf
(gsub("[A-Z|a-z]|[.]","",HGVSp_Short) %in%
# if overlaps the hotspot Amino_Acid_Position
hotspot_database_2017_tp53_snv$Amino_Acid_Position) ~ 1,
TRUE ~ 0),
activating = case_when(
# strip REF and Variant AA to get Amino_Acid_Position in consensus maf
(gsub("[A-Z|a-z]|[.]","",HGVSp_Short) %in%
# if overlaps the activating Amino_Acid_Position
c("273","248")) ~ 1,
TRUE ~ 0
)
)
Gather SNV,CNV, classifier and cancer_predisposition values
We will gather values per sample_id since DNA and RNA samples can be only matched by sample_id
tp53_alterations_dna <- histology %>%
filter(experimental_strategy != "RNA-Seq") %>%
# merge cell line composition
left_join(cell_line_composition) %>%
# filter for Tumors
filter(sample_type=="Tumor") %>%
# join consensus calls overlapping hotspots
left_join(consensus_tp53_snv_indel,
by=c("Kids_First_Biospecimen_ID"="Tumor_Sample_Barcode")) %>%
# join filtered cnv losses
left_join(cnv_domain_overlap,by=c("Kids_First_Biospecimen_ID"="biospecimen_id",
"sample_id")) %>%
# join SV
left_join(sv_overlap,by=c("Kids_First_Biospecimen_ID"="Kids.First.Biospecimen.ID.Tumor",
"sample_id")) %>%
dplyr::rename(Kids_First_Biospecimen_ID_DNA = Kids_First_Biospecimen_ID)
Joining, by = c("Kids_First_Biospecimen_ID", "sample_id", "Kids_First_Participant_ID", "cell_line_composition", "aliquot_id")
tp53_alterations_score <- tp53_alterations_dna %>%
# add classifier score
full_join(classifier_score,by=c("sample_id",
"cell_line_composition",
"cancer_predispositions",
"sample_type",
"Kids_First_Participant_ID")) %>%
# select useful columns
select(sample_id,Kids_First_Biospecimen_ID_RNA,Kids_First_Biospecimen_ID_DNA,
tp53_score,cancer_predispositions,HGVSp_Short,copy_number,SV.type,hotspot,activating) %>%
arrange(sample_id) %>%
unique() %>%
replace_na(list(hotspot = 0, activating = 0))
tp53_alterations_score
Convert the above to wide format
tp53_alterations_score_wide <- tp53_alterations_score %>%
# group
group_by(sample_id,
Kids_First_Biospecimen_ID_DNA,
Kids_First_Biospecimen_ID_RNA,
cancer_predispositions,
tp53_score) %>%
#
# sample_id as rows add SNV counts,CNV loss counts
# and hotspot/activating mutation annotation
summarise(
# summarize length of SNV and CNV alterations
SNV_indel_counts= length(unique(HGVSp_Short[!is.na(HGVSp_Short)])),
CNV_loss_counts = length(unique(copy_number[!is.na(copy_number)])),
SV_counts = length(unique(SV.type[!is.na(SV.type)])),
HGVSp_Short = toString(unique(HGVSp_Short)),
CNV_loss_evidence = toString(unique(copy_number)),
SV_type = toString(unique(SV.type)),
# summarize unique hotspot values per sample_id
hotspot = max(unique(hotspot[!is.na(hotspot) ])),
activating = max(unique(activating[!is.na(activating)])))
tp53_alterations_score_wide
Add annotation
As discussed above we want to annotate TP53 mutants that are putative loss-of-function OR gain-of-function.
tp53_alterations_score_wide <- tp53_alterations_score_wide %>%
mutate(
# add tp53_altered annotation column
tp53_altered =
case_when(
# when activating == 0
activating == 0 &
# check if mutated variant AA position overlaps hotspot database
( hotspot == 1 |
# check if sample_id has SNV+(CNV|SV) mutation
# suggesting both alleles are mutated
(SNV_indel_counts >= 1 &
(CNV_loss_counts >=1 | SV_counts >=1) ) |
# check if more than 1 SNV is present
# suggesting both alleles are mutated
SNV_indel_counts > 1 |
# check if SNV mutant and has
# Li-Fraumeni syndrome suggesting
# germline TP53 mutant as cancer predisposition
(SNV_indel_counts >= 1 &
grepl("Li-Fraumeni syndrome",cancer_predispositions)) |
# check if CNV|SV mutant and has
# Li-Fraumeni syndrome suggesting
# germline TP53 mutant as cancer predisposition
((CNV_loss_counts >= 1 | SV_counts >=1 ) &
grepl("Li-Fraumeni syndrome",cancer_predispositions)) |
# check if tp53 inactivating score for RNA_Seq is greater that 0.5
# and has Li-Fraumeni syndrome suggesting
# germline TP53 mutant as cancer predisposition
(tp53_score > 0.5 &
grepl("Li-Fraumeni syndrome",cancer_predispositions)) |
# check if tp53 inactivating score for RNA_Seq is greater that 0.5
# and has 1 or more SNV| (CNV|SV) loss in TP53
(tp53_score > 0.5 &
(SNV_indel_counts >= 1 | (CNV_loss_counts >=1 | SV_counts >=1 )))
) ~ "loss",
# when activating == 1
activating == 1 ~ "activated",
# if no evidence supports TP53 inativation or activation
TRUE ~ "Other"
)
)
Explore distribution of tp53_altered status vs tp53 inactivation scores
ggplot(tp53_alterations_score_wide, aes(x = factor(tp53_altered), y = tp53_score)) +
geom_violin()+
geom_jitter(alpha = 0.5, width = 0.2) +
stat_compare_means() +
theme_bw() +
ggtitle("Distribution of scores across tp53 altered status") +
theme(axis.text.x = element_text(angle = 60, hjust = 1))+
xlab("tp53 altered status")
Warning: Removed 96 rows containing non-finite values (stat_ydensity).
Warning: Removed 96 rows containing non-finite values (stat_compare_means).
Warning: Removed 96 rows containing missing values (geom_point).

TP53 mutants with gain-of-function (activating) mutants promote tumorigenesis according to literature. Reference and reference have similar distribution of classifier scores to that of potential loss-of-function (multi-allelic) Tp53 mutations.
“Other” annotated samples either don’t have any high-confidence loss/gain TP53 SNVs nor CNV losses OR DNA sample is not available.
Expression profile for activating vs loss TP53 status
expr_combined <- readRDS(file.path(data_dir,"gene-counts-rsem-expected_count-collapsed.rds"))
(stranded) Expression profile for activating vs loss TP53 status
# subset to stranded samples
stranded_samples <- clinical %>%
filter(sample_type == "Tumor",
experimental_strategy == "RNA-Seq",
cohort %in% c("PBTA","GMKF","TARGET"),
RNA_library == "stranded") %>%
pull(Kids_First_Biospecimen_ID) %>% unique()
stranded <- expr_combined %>% select(stranded_samples)
# subset to TP53
subset_stranded <- t(stranded)[,"TP53"]
stranded_tp53 <- as.data.frame(subset_stranded) %>%
dplyr::rename(TP53=subset_stranded) %>%
rownames_to_column() %>%
left_join(tp53_alterations_score_wide,by=c("rowname"="Kids_First_Biospecimen_ID_RNA")) %>%
# remove polya from matched scores
filter(tp53_altered %in% c("activated","loss"))
Plot distribution of TP53 gene expression (stranded)
ggplot(stranded_tp53, aes(x = factor(tp53_altered), y = TP53)) +
geom_violin()+
geom_jitter(alpha = 0.5, width = 0.2) +
stat_compare_means() +
theme_bw() +
ggtitle("Distribution of TP53 expression across tp53 altered status") +
theme(axis.text.x = element_text(angle = 60, hjust = 1))+
xlab("tp53 altered status")

(poly-a) Expression profile for activating vs loss TP53 status
# subset to polya samples
polya_samples <- clinical %>%
filter(sample_type == "Tumor",
experimental_strategy == "RNA-Seq",
cohort %in% c("PBTA","GMKF","TARGET"),
RNA_library == "poly-A") %>%
pull(Kids_First_Biospecimen_ID) %>% unique()
polya <- expr_combined %>% select(polya_samples)
# subset to TP53
subset_polya<- t(polya)[,"TP53"]
polya_tp53<- as.data.frame(subset_polya) %>%
dplyr::rename(TP53=subset_polya) %>%
rownames_to_column() %>%
left_join(tp53_alterations_score_wide,by=c("rowname"="Kids_First_Biospecimen_ID_RNA")) %>%
filter(tp53_altered %in% c("activated","loss"))
Plot distribution of TP53 gene expression (poly-a)
ggplot(polya_tp53, aes(x = factor(tp53_altered), y = TP53)) +
geom_violin()+
geom_jitter(alpha = 0.5, width = 0.2) +
stat_compare_means() +
theme_bw() +
ggtitle("Distribution of TP53 expression across tp53 altered status") +
theme(axis.text.x = element_text(angle = 60, hjust = 1))+
xlab("tp53 altered status")

(poly-a stranded) Expression profile for activating vs loss TP53 status
# subset to polya samples
polya_stranded_samples <- clinical %>%
filter(sample_type == "Tumor",
experimental_strategy == "RNA-Seq",
cohort %in% c("PBTA","GMKF","TARGET"),
RNA_library == "poly-A stranded") %>%
pull(Kids_First_Biospecimen_ID) %>% unique()
polya_stranded <- expr_combined %>% select(polya_stranded_samples)
# subset to TP53
subset_polya_stranded<- t(polya_stranded)[,"TP53"]
polya_stranded_tp53<- as.data.frame(subset_polya_stranded) %>%
dplyr::rename(TP53=subset_polya_stranded) %>%
rownames_to_column() %>%
left_join(tp53_alterations_score_wide,by=c("rowname"="Kids_First_Biospecimen_ID_RNA")) %>%
filter(tp53_altered %in% c("activated","loss"))
Plot distribution of TP53 gene expression (poly-a)
ggplot(polya_stranded_tp53, aes(x = factor(tp53_altered), y = TP53)) +
geom_violin()+
geom_jitter(alpha = 0.5, width = 0.2) +
stat_compare_means() +
theme_bw() +
ggtitle("Distribution of TP53 expression across tp53 altered status") +
theme(axis.text.x = element_text(angle = 60, hjust = 1))+
xlab("tp53 altered status")

Check if other cancer predisposition have high TP53 inactivation scores
ggplot(tp53_alterations_score_wide, aes(x = factor(tp53_altered), y = tp53_score)) +
geom_violin()+
geom_jitter(alpha = 0.5, width = 0.2) +
theme_bw() +
ggtitle("Distribution of scores across tp53 altered status") +
theme(axis.text.x = element_text(angle = 60, hjust = 1))+
xlab("tp53 altered status") +
facet_wrap(.~cancer_predispositions)
Warning: Removed 96 rows containing non-finite values (stat_ydensity).
Warning in max(data$density): no non-missing arguments to max; returning -
Inf
Warning in max(data$density): no non-missing arguments to max; returning -
Inf
Warning: Removed 96 rows containing missing values (geom_point).

Some NF-1 and/or Other inherited conditions NOS have high scoring TP53 mutants as well.
Interesting 2 bs ids annotated with Li-Fraumeni synfrome pre-disposition have very low tp53 classifier scores.
tp53_alterations_score_wide %>%
filter(cancer_predispositions == "Li-Fraumeni syndrome",
tp53_score < 0.5)
Save file
Adding DNA and RNA biospecimen ids to tp53_alterations_score_wide and saving
tp53_alterations_score_wide %>%
write_tsv(file.path(results_dir,"tp53_altered_status.tsv"))
LS0tCnRpdGxlOiAiVHA1MyBTTlYgaG90c3BvdHMiCmF1dGhvcjogIksgUyBHYW9ua2FyIChEM0IpIgpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sKcGFyYW1zOgogIGJhc2VfcnVuOgogICAgbGFiZWw6ICIxLzAgdG8gcnVuIHdpdGggYmFzZSBoaXN0b2xvZ3kiCiAgICB2YWx1ZTogMAogICAgaW5wdXQ6IGludGVnZXIKLS0tCgpJbiB0aGlzIG5vdGVib29rIHdlIHdpbGwgYWRkIFRQNTMgYWx0ZXJhdGlvbiBzdGF0dXMgYXMgZGlzY3Vzc2VkIGluIFsjODM3XShodHRwczovL2dpdGh1Yi5jb20vQWxleHNMZW1vbmFkZS9PcGVuUEJUQS1hbmFseXNpcy9pc3N1ZXMvODM3KToKCioqVFA1MyBhbHRlcmVkIC0gbG9zcyoqLCBpZiA6CgogLSBBIHNhbXBsZSBjb250YWlucyBhIFRQNTMgaG90c3BvdCBTTlYgbXV0YXRpb24uIChDYW5jZXIgaG90c3BvdCBkYXRhYmFzZSBhbmQgZG93bmxvYWRhYmxlIGZpbGUgYXZhaWxhYmxlKS4gUGxlYXNlIGFsc28gY3Jvc3NjaGVjayB0aGF0IGFsbCBtdXRhdGlvbnMgZnJvbSB0aGlzIHRhYmxlIGFyZSBpbmNsdWRlZC4KIC0gQSBzYW1wbGUgY29udGFpbnMgdHdvIFRQNTMgYWx0ZXJhdGlvbnMsIHN1Z2dlc3RpbmcgKGJ1dCBub3QgY29uZmlybWluZykgdGhhdCBib3RoIGFsbGVsZXMgYXJlIGFmZmVjdGVkIChTTlYrU05WLCAoQ05WIG9yIFNWKSArIFNOVikuCiAtIEEgc2FtcGxlIGNvbnRhaW5zIG9uZSBhbHRlcmF0aW9uIChTTlYgb3IgKENOViBvciBTVikpICsgaGFzIGNhbmNlcl9wcmVkaXNwb3NpdGlvbnMgPT0gIkxpLUZyYXVtZW5pIHN5bmRyb21lIiwgc3VnZ2VzdGluZyB0aGVyZSBpcyBhIGdlcm1saW5lIHZhcmlhbnQgaW4gYWRkaXRpb24gdG8gdGhlIHNvbWF0aWMgdmFyaWFudCB3ZSBvYnNlcnZlLgogLSBBIHNhbXBsZSBkb2VzIG5vdCBoYXZlIGEgVFA1MyBhbHRlcmF0aW9ucywgYnV0IGhhcyBjYW5jZXJfcHJlZGlzcG9zaXRpb25zID09ICJMaS1GcmF1bWVuaSBzeW5kcm9tZSIgYW5kIFRQNTMgY2xhc3NpZmllciBzY29yZSBmb3IgbWF0Y2hlZCBSTkEtU2VxID4gMC41IChvciBoaWdoZXIgY3V0b2ZmIHdlIGRlY2lkZSB1cG9uIGxhdGVyKS4KIAogTm90ZTogQ05WIGFuZCBTViB3aWxsIGJlIGNvbnNpZGVyZWQgYXMgdGhlIHNhbWUgZXZlbnQsIGJ1dCB3ZSB3aWxsIGJlIGFkZGluZyBTVl9jb3VudHMgYW5kIFNWLnR5cGUgdG8gdGhlIGFsdGVyZWQgc3RhdHVzIG91dHB1dCBmaWxlCiAKKipUUDUzIGFsdGVyZWQgLSBhY3RpdmF0ZWQqKiwgaWY6CgotIEEgc2FtcGxlIGNvbnRhaW5zIG9uZSBvZiB0aGUgdHdvIFRQNTMgYWN0aXZhdGluZyBtdXRhdGlvbnMgcC5SMjczQyBhbmQgcC5SMjQ4Vy4gW0Bkb2k6MTAuMTAzOC9uZzA1OTMtNDJdCgojIyBTZXR1cApgYGB7cn0KbGlicmFyeSgiZ2dwdWJyIikKbGlicmFyeSgiZ2d0aGVtZXMiKQpsaWJyYXJ5KCJ0aWR5dmVyc2UiKQpsaWJyYXJ5KCJnZ2ZvcnRpZnkiKQpsaWJyYXJ5KCJicm9vbSIpCgojIHJvb3RkaXIKcm9vdF9kaXIgPC0gcnByb2pyb290OjpmaW5kX3Jvb3QocnByb2pyb290OjpoYXNfZGlyKCIuZ2l0IikpCmRhdGFfZGlyIDwtIGZpbGUucGF0aChyb290X2RpciwgImRhdGEiKQoKaW5wdXRfZGlyIDwtIGZpbGUucGF0aChyb290X2RpciwKICAgICAgICAgICAgICAgICAgICAgICAgICJhbmFseXNlcyIsCiAgICAgICAgICAgICAgICAgICAgICAgICAidHA1M19uZjFfc2NvcmUiLAogICAgICAgICAgICAgICAgICAgICAgICAgImlucHV0IikKIyBjZWxsIGNvbXBvc2l0aW9uCmNlbGxfbGluZV9jb21wb3NpdGlvbiA8LSByZWFkX3RzdihmaWxlLnBhdGgoIi4uIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAibW9sZWN1bGFyLXN1YnR5cGluZy1IR0ciLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJpbnB1dCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgImNlbGwtbGluZS1jb21wb3NpdGlvbi50c3YiKSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGNvbF90eXBlcyA9IAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICByZWFkcjo6Y29scyggYWxpcXVvdF9pZCAgPSByZWFkcjo6Y29sX2NoYXJhY3RlcigpKSkgCgojIGhpc3RvbG9neQppZiAoIHBhcmFtcyRiYXNlX3J1biA9PTAgKXsKICBjbGluaWNhbDwtcmVhZC5kZWxpbShmaWxlLnBhdGgoZGF0YV9kaXIsImhpc3RvbG9naWVzLnRzdiIpLCBzdHJpbmdzQXNGYWN0b3JzID0gRkFMU0UpCn0gZWxzZXsKICBjbGluaWNhbDwtcmVhZC5kZWxpbShmaWxlLnBhdGgoZGF0YV9kaXIsImhpc3RvbG9naWVzLWJhc2UudHN2IiksIHN0cmluZ3NBc0ZhY3RvcnMgPSBGQUxTRSkgIAp9CgpoaXN0b2xvZ3kgPC0gY2xpbmljYWwgJT4lCiAgIHNlbGVjdCgiS2lkc19GaXJzdF9CaW9zcGVjaW1lbl9JRCIsCiAgICAgICAgICAic2FtcGxlX2lkIiwKICAgICAgICAgICJjYW5jZXJfcHJlZGlzcG9zaXRpb25zIiwKICAgICAgICAgICJzYW1wbGVfdHlwZSIsCiAgICAgICAgICAiZXhwZXJpbWVudGFsX3N0cmF0ZWd5IiwKICAgICAgICAgICJjb21wb3NpdGlvbiIpICU+JQogICAjIG1lcmdlIGNlbGwgbGluZSBjb21wb3NpdGlvbgogIGxlZnRfam9pbihjZWxsX2xpbmVfY29tcG9zaXRpb24pCiAKCiMgQ2FuY2VyIGRhdGFiYXNlCmhvdHNwb3RfZGF0YWJhc2VfMjAxN190cDUzX3NudiA8LSByZWFkeGw6OnJlYWRfeGxzKGZpbGUucGF0aChpbnB1dF9kaXIsImhvdHNwb3RzX3YyLnhscyIpLHNoZWV0ID0gMSkgJT4lCiAgZmlsdGVyKEh1Z29fU3ltYm9sID09ICJUUDUzIikKaG90c3BvdF9kYXRhYmFzZV8yMDE3X3RwNTNfaW5kZWwgPC0gcmVhZHhsOjpyZWFkX3hscyhmaWxlLnBhdGgoaW5wdXRfZGlyLCJob3RzcG90c192Mi54bHMiKSxzaGVldCA9IDIpICU+JQogIGZpbHRlcihIdWdvX1N5bWJvbCA9PSAiVFA1MyIpCgojIHA1MywgcDYzIGFuZCBwNzMgZnVuY3Rpb25hbCBzaXRlcwpmdW5jdGlvbmFsX3NpdGVzX3RwNTMgPC0gcmVhZF90c3YoZmlsZS5wYXRoKGlucHV0X2RpciwiaG90c3BvdF9DaGVuXzIwMDYudHN2IikpIAoKcmVzdWx0c19kaXIgPC0gZmlsZS5wYXRoKHJvb3RfZGlyLAogICAgICAgICAgICAgICAgICAgICAgICAgImFuYWx5c2VzIiwKICAgICAgICAgICAgICAgICAgICAgICAgICJ0cDUzX25mMV9zY29yZSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAicmVzdWx0cyIpCgppZiAoIWRpci5leGlzdHMocmVzdWx0c19kaXIpKSB7CiAgZGlyLmNyZWF0ZShyZXN1bHRzX2RpcikKfQoKIyBjbGFzc2lmaWVyIHNjb3JlcwpjbGFzc2lmaWVyX3Njb3JlIDwtIHJlYWRfdHN2KAogIGZpbGUucGF0aCgKICAgIHJlc3VsdHNfZGlyLAogICAgImdlbmUtZXhwcmVzc2lvbi1yc2VtLXRwbS1jb2xsYXBzZWRfY2xhc3NpZmllcl9zY29yZXMudHN2IikpJT4lCiAgc2VsZWN0KHNhbXBsZV9pZCx0cDUzX3Njb3JlKSAlPiUKICBkcGx5cjo6cmVuYW1lKEtpZHNfRmlyc3RfQmlvc3BlY2ltZW5fSURfUk5BPXNhbXBsZV9pZCkgJT4lCiAgbGVmdF9qb2luKGhpc3RvbG9neSAsYnk9YygiS2lkc19GaXJzdF9CaW9zcGVjaW1lbl9JRF9STkEiPSJLaWRzX0ZpcnN0X0Jpb3NwZWNpbWVuX0lEIikpIAoKIyBDb3B5IG51bWJlciBwZXIgc2FtcGxlIAojIG92ZXJsYXBwaW5nIGZ1bmN0aW9uYWwgZG9tYWluCmNudl9kb21haW5fb3ZlcmxhcCA8LSByZWFkX3RzdigKICBmaWxlLnBhdGgoCiAgICByZXN1bHRzX2RpciwKICAgICJsb3NzX292ZXJsYXBfZG9tYWluc190cDUzLnRzdiIpKQoKIyBTdHJ1Y3R1cmFsIHZhcmlhbnQgb3ZlcmxhcHBpbmcgCiMgb3Igd2l0aGluIGdlbmUgbG9jdXMgb2YgVFA1Mwpzdl9vdmVybGFwIDwtIHJlYWRfdHN2KAogIGZpbGUucGF0aCgKICAgIHJlc3VsdHNfZGlyLAogICAgInN2X292ZXJsYXBfdHA1My50c3YiKSkKCmBgYAoKIyMjIENoZWNrIGlmIGFsbCBmdW5jdGlvbiBzaXRlcyBhcmUgaW4gaG90c3BvdHMKCkFsbCBmdW5jdGlvbmFsIHNpdGVzIGFyZSBqdXN0IDEgYmFzZSBzbyBjaGVja2luZyBpbiBgaG90c3BvdF9kYXRhYmFzZV8yMDE3X3NudmAKCmBgYHtyfQoKZnVuY3Rpb25hbF9zaXRlc190cDUzICU+JQogIGZpbHRlcighZ3N1YigiW0EtWnxhLXpdIiwiIixwNTMpICVpbiUgaG90c3BvdF9kYXRhYmFzZV8yMDE3X3RwNTNfc252JEFtaW5vX0FjaWRfUG9zaXRpb24gKQoKYGBgCjIgZnVuY3Rpb25hbCBzaXRlcyBhcmUgbWlzc2luZyBpbiBob3RzcG90cyBkYXRhYmFzZS4KCiMjIFNOViBmb3IgVFA1MwoKUmVtb3ZpbmcgU2lsZW50IG9yIEludHJvbiBjbGFzc2lmaWVkIHZhcmlhbnRzIHRvIGNhcHR1cmUgb25seSBwdXRhdGl2ZSBkYW1hZ2luZyBtdXRhdGlvbnMKCmBgYHtyfQpjb25zZW5zdXNfdHA1M19zbnZfaW5kZWwgPC0gZGF0YS50YWJsZTo6ZnJlYWQoCiAgZmlsZS5wYXRoKGRhdGFfZGlyLCJzbnYtY29uc2Vuc3VzLXBsdXMtaG90c3BvdHMubWFmLnRzdi5neiIpLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHNlbGVjdCA9IGMoIkNocm9tb3NvbWUiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlN0YXJ0X1Bvc2l0aW9uIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJFbmRfUG9zaXRpb24iLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlN0cmFuZCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiVmFyaWFudF9DbGFzc2lmaWNhdGlvbiIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiVHVtb3JfU2FtcGxlX0JhcmNvZGUiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkh1Z29fU3ltYm9sIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJIR1ZTcF9TaG9ydCIpLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGRhdGEudGFibGUgPSBGQUxTRSkgJT4lCiAgZmlsdGVyKEh1Z29fU3ltYm9sID09ICJUUDUzIikgJT4lCiAgZmlsdGVyKCEoVmFyaWFudF9DbGFzc2lmaWNhdGlvbiAlaW4lIGMoIlNpbGVudCIsICJJbnRyb24iLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICMgcmVtb3ZlIG90aGVyIG5vbi1hbWlubyBhY2lkIFNOVnMKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiMydGbGFuayIgLCI1J0ZsYW5rIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiMydVVFIiLCAiNSdVVFIiICApKSkKYGBgCgojIyMgIEdhdGhlciBhbm5vdGF0aW9uIGZvciBTTlYgCgpob3RzcG90cyA6IG92ZXJsYXBzIEFBIHBvc2l0aW9uIHdoaWNoIGFyZSBzdGF0aXN0aWNhbGx5IHNpZ25pZmljYW50IFNOVi9JbmRlbHMKYWN0aXZhdGluZyA6IG92ZXJsYXBzIEFBIHBvc2l0aW9uIHdoaWNoIGFyZSBmb3VuZCB0byBhY3QgYXMgZ2Fpbi1vZi1mdW5jdGlvbiBhY2NvcmRpbmcgdG8gbGl0ZXJhdHVyZQoKYGBge3J9Cgpjb25zZW5zdXNfdHA1M19zbnZfaW5kZWwgPC0gY29uc2Vuc3VzX3RwNTNfc252X2luZGVsICU+JQogIG11dGF0ZSgKICAgIGhvdHNwb3QgPSBjYXNlX3doZW4oCiAgICAgICMgc3RyaXAgUkVGIGFuZCBWYXJpYW50IEFBIHRvIGdldCBBbWlub19BY2lkX1Bvc2l0aW9uIGluIGNvbnNlbnN1cyBtYWYKICAgICAgKGdzdWIoIltBLVp8YS16XXxbLl0iLCIiLEhHVlNwX1Nob3J0KSAlaW4lIAogICAgICAgICAjIGlmIG92ZXJsYXBzIHRoZSBob3RzcG90IEFtaW5vX0FjaWRfUG9zaXRpb24gCiAgICAgICAgIGhvdHNwb3RfZGF0YWJhc2VfMjAxN190cDUzX3NudiRBbWlub19BY2lkX1Bvc2l0aW9uKSAgfiAxLAogICAgICBUUlVFIH4gMCksCiAgICBhY3RpdmF0aW5nID0gY2FzZV93aGVuKAogICAgICAjIHN0cmlwIFJFRiBhbmQgVmFyaWFudCBBQSB0byBnZXQgQW1pbm9fQWNpZF9Qb3NpdGlvbiBpbiBjb25zZW5zdXMgbWFmCiAgICAgIChnc3ViKCJbQS1afGEtel18Wy5dIiwiIixIR1ZTcF9TaG9ydCkgJWluJSAKICAgICAgICAgIyBpZiBvdmVybGFwcyB0aGUgYWN0aXZhdGluZyAgQW1pbm9fQWNpZF9Qb3NpdGlvbiAKICAgICAgICAgYygiMjczIiwiMjQ4IikpICB+IDEsCiAgICAgIFRSVUUgfiAwCiAgICApCiAgKQoKYGBgCgojIyBHYXRoZXIgU05WLENOViwgY2xhc3NpZmllciBhbmQgY2FuY2VyX3ByZWRpc3Bvc2l0aW9uIHZhbHVlcwoKV2Ugd2lsbCBnYXRoZXIgdmFsdWVzIHBlciBzYW1wbGVfaWQgc2luY2UgRE5BIGFuZCBSTkEgc2FtcGxlcyBjYW4gYmUgb25seSBtYXRjaGVkIGJ5IGBzYW1wbGVfaWRgCgoKYGBge3J9Cgp0cDUzX2FsdGVyYXRpb25zX2RuYSA8LSBoaXN0b2xvZ3kgJT4lCiAgZmlsdGVyKGV4cGVyaW1lbnRhbF9zdHJhdGVneSAhPSAiUk5BLVNlcSIpICU+JQogICMgbWVyZ2UgY2VsbCBsaW5lIGNvbXBvc2l0aW9uCiAgbGVmdF9qb2luKGNlbGxfbGluZV9jb21wb3NpdGlvbikgJT4lCiAgIyBmaWx0ZXIgZm9yIFR1bW9ycwogIGZpbHRlcihzYW1wbGVfdHlwZT09IlR1bW9yIikgJT4lCiAgIyBqb2luIGNvbnNlbnN1cyBjYWxscyBvdmVybGFwcGluZyBob3RzcG90cwogIGxlZnRfam9pbihjb25zZW5zdXNfdHA1M19zbnZfaW5kZWwsIAogICAgICAgICAgICBieT1jKCJLaWRzX0ZpcnN0X0Jpb3NwZWNpbWVuX0lEIj0iVHVtb3JfU2FtcGxlX0JhcmNvZGUiKSkgJT4lCiAgIyBqb2luIGZpbHRlcmVkIGNudiBsb3NzZXMKICBsZWZ0X2pvaW4oY252X2RvbWFpbl9vdmVybGFwLGJ5PWMoIktpZHNfRmlyc3RfQmlvc3BlY2ltZW5fSUQiPSJiaW9zcGVjaW1lbl9pZCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJzYW1wbGVfaWQiKSkgJT4lCiAgIyBqb2luIFNWIAogIGxlZnRfam9pbihzdl9vdmVybGFwLGJ5PWMoIktpZHNfRmlyc3RfQmlvc3BlY2ltZW5fSUQiPSJLaWRzLkZpcnN0LkJpb3NwZWNpbWVuLklELlR1bW9yIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICJzYW1wbGVfaWQiKSkgJT4lCiAgZHBseXI6OnJlbmFtZShLaWRzX0ZpcnN0X0Jpb3NwZWNpbWVuX0lEX0ROQSA9IEtpZHNfRmlyc3RfQmlvc3BlY2ltZW5fSUQpIAoKdHA1M19hbHRlcmF0aW9uc19zY29yZSA8LSB0cDUzX2FsdGVyYXRpb25zX2RuYSAlPiUKICAjIGFkZCBjbGFzc2lmaWVyIHNjb3JlCiAgZnVsbF9qb2luKGNsYXNzaWZpZXJfc2NvcmUsYnk9Yygic2FtcGxlX2lkIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJjZWxsX2xpbmVfY29tcG9zaXRpb24iLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgImNhbmNlcl9wcmVkaXNwb3NpdGlvbnMiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgInNhbXBsZV90eXBlIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJLaWRzX0ZpcnN0X1BhcnRpY2lwYW50X0lEIikpICU+JQogICMgc2VsZWN0IHVzZWZ1bCBjb2x1bW5zCiAgc2VsZWN0KHNhbXBsZV9pZCxLaWRzX0ZpcnN0X0Jpb3NwZWNpbWVuX0lEX1JOQSxLaWRzX0ZpcnN0X0Jpb3NwZWNpbWVuX0lEX0ROQSwKICAgICAgICAgdHA1M19zY29yZSxjYW5jZXJfcHJlZGlzcG9zaXRpb25zLEhHVlNwX1Nob3J0LGNvcHlfbnVtYmVyLFNWLnR5cGUsaG90c3BvdCxhY3RpdmF0aW5nKSAlPiUKICBhcnJhbmdlKHNhbXBsZV9pZCkgJT4lCiAgdW5pcXVlKCkgJT4lCiAgcmVwbGFjZV9uYShsaXN0KGhvdHNwb3QgPSAwLCBhY3RpdmF0aW5nID0gMCkpCgp0cDUzX2FsdGVyYXRpb25zX3Njb3JlIAoKYGBgCgpDb252ZXJ0IHRoZSBhYm92ZSB0byB3aWRlIGZvcm1hdAoKYGBge3J9Cgp0cDUzX2FsdGVyYXRpb25zX3Njb3JlX3dpZGUgPC0gdHA1M19hbHRlcmF0aW9uc19zY29yZSAlPiUKICAjIGdyb3VwIAogIGdyb3VwX2J5KHNhbXBsZV9pZCwKICAgICAgICAgICBLaWRzX0ZpcnN0X0Jpb3NwZWNpbWVuX0lEX0ROQSwKICAgICAgICAgICBLaWRzX0ZpcnN0X0Jpb3NwZWNpbWVuX0lEX1JOQSwKICAgICAgICAgICBjYW5jZXJfcHJlZGlzcG9zaXRpb25zLAogICAgICAgICAgIHRwNTNfc2NvcmUpICU+JQogICMgCiAgIyBzYW1wbGVfaWQgYXMgcm93cyBhZGQgU05WIGNvdW50cyxDTlYgbG9zcyBjb3VudHMgCiAgIyBhbmQgaG90c3BvdC9hY3RpdmF0aW5nIG11dGF0aW9uIGFubm90YXRpb24gIAogIHN1bW1hcmlzZSgKICAgICMgc3VtbWFyaXplIGxlbmd0aCBvZiBTTlYgYW5kIENOViBhbHRlcmF0aW9ucyAKICAgIFNOVl9pbmRlbF9jb3VudHM9IGxlbmd0aCh1bmlxdWUoSEdWU3BfU2hvcnRbIWlzLm5hKEhHVlNwX1Nob3J0KV0pKSwgCiAgICBDTlZfbG9zc19jb3VudHMgPSBsZW5ndGgodW5pcXVlKGNvcHlfbnVtYmVyWyFpcy5uYShjb3B5X251bWJlcildKSksCiAgICBTVl9jb3VudHMgPSBsZW5ndGgodW5pcXVlKFNWLnR5cGVbIWlzLm5hKFNWLnR5cGUpXSkpLAogICAgSEdWU3BfU2hvcnQgPSB0b1N0cmluZyh1bmlxdWUoSEdWU3BfU2hvcnQpKSwKICAgIENOVl9sb3NzX2V2aWRlbmNlID0gdG9TdHJpbmcodW5pcXVlKGNvcHlfbnVtYmVyKSksCiAgICBTVl90eXBlID0gdG9TdHJpbmcodW5pcXVlKFNWLnR5cGUpKSwKICAgICMgc3VtbWFyaXplIHVuaXF1ZSBob3RzcG90IHZhbHVlcyBwZXIgc2FtcGxlX2lkCiAgICBob3RzcG90ID0gbWF4KHVuaXF1ZShob3RzcG90WyFpcy5uYShob3RzcG90KSBdKSksCiAgICBhY3RpdmF0aW5nID0gbWF4KHVuaXF1ZShhY3RpdmF0aW5nWyFpcy5uYShhY3RpdmF0aW5nKV0pKSkgCgp0cDUzX2FsdGVyYXRpb25zX3Njb3JlX3dpZGUKCmBgYAoKIyMjIEFkZCBhbm5vdGF0aW9uCgpBcyBkaXNjdXNzZWQgYWJvdmUgd2Ugd2FudCB0byBhbm5vdGF0ZSBUUDUzIG11dGFudHMgdGhhdCBhcmUgcHV0YXRpdmUgbG9zcy1vZi1mdW5jdGlvbiBPUiAgZ2Fpbi1vZi1mdW5jdGlvbi4KCmBgYHtyfQoKdHA1M19hbHRlcmF0aW9uc19zY29yZV93aWRlIDwtIHRwNTNfYWx0ZXJhdGlvbnNfc2NvcmVfd2lkZSAlPiUKICBtdXRhdGUoCiAgICAjIGFkZCB0cDUzX2FsdGVyZWQgYW5ub3RhdGlvbiBjb2x1bW4KICAgIHRwNTNfYWx0ZXJlZCA9IAogICAgICBjYXNlX3doZW4oCiAgICAgICAgIyB3aGVuIGFjdGl2YXRpbmcgPT0gMAogICAgICAgIGFjdGl2YXRpbmcgPT0gMCAmCiAgICAgICAgICAjIGNoZWNrIGlmIG11dGF0ZWQgdmFyaWFudCBBQSBwb3NpdGlvbiBvdmVybGFwcyBob3RzcG90IGRhdGFiYXNlIAogICAgICAgICAgKCBob3RzcG90ID09IDEgIHwKICAgICAgICAgICAgICAjIGNoZWNrIGlmIHNhbXBsZV9pZCBoYXMgU05WKyhDTlZ8U1YpIG11dGF0aW9uCiAgICAgICAgICAgICAgIyBzdWdnZXN0aW5nIGJvdGggYWxsZWxlcyBhcmUgbXV0YXRlZAogICAgICAgICAgICAgIChTTlZfaW5kZWxfY291bnRzID49IDEgJiAKICAgICAgICAgICAgICAgICAoQ05WX2xvc3NfY291bnRzID49MSB8IFNWX2NvdW50cyA+PTEpICkgfAogICAgICAgICAgICAgICMgY2hlY2sgaWYgbW9yZSB0aGFuIDEgU05WIGlzIHByZXNlbnQKICAgICAgICAgICAgICAjIHN1Z2dlc3RpbmcgYm90aCBhbGxlbGVzIGFyZSBtdXRhdGVkCiAgICAgICAgICAgICAgU05WX2luZGVsX2NvdW50cyA+IDEgfAogICAgICAgICAgICAgICMgY2hlY2sgaWYgU05WIG11dGFudCBhbmQgaGFzCiAgICAgICAgICAgICAgIyBMaS1GcmF1bWVuaSBzeW5kcm9tZSBzdWdnZXN0aW5nCiAgICAgICAgICAgICAgIyBnZXJtbGluZSBUUDUzIG11dGFudCBhcyBjYW5jZXIgcHJlZGlzcG9zaXRpb24gCiAgICAgICAgICAgICAgKFNOVl9pbmRlbF9jb3VudHMgPj0gMSAmIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgIGdyZXBsKCJMaS1GcmF1bWVuaSBzeW5kcm9tZSIsY2FuY2VyX3ByZWRpc3Bvc2l0aW9ucykpIHwKICAgICAgICAgICAgICAjIGNoZWNrIGlmIENOVnxTViBtdXRhbnQgIGFuZCBoYXMKICAgICAgICAgICAgICAjIExpLUZyYXVtZW5pIHN5bmRyb21lIHN1Z2dlc3RpbmcKICAgICAgICAgICAgICAjIGdlcm1saW5lIFRQNTMgbXV0YW50IGFzIGNhbmNlciBwcmVkaXNwb3NpdGlvbiAKICAgICAgICAgICAgICAoKENOVl9sb3NzX2NvdW50cyA+PSAxIHwgU1ZfY291bnRzID49MSApICYgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgZ3JlcGwoIkxpLUZyYXVtZW5pIHN5bmRyb21lIixjYW5jZXJfcHJlZGlzcG9zaXRpb25zKSkgfAogICAgICAgICAgICAgICMgY2hlY2sgaWYgdHA1MyBpbmFjdGl2YXRpbmcgc2NvcmUgZm9yIFJOQV9TZXEgaXMgZ3JlYXRlciB0aGF0IDAuNQogICAgICAgICAgICAgICMgYW5kIGhhcyBMaS1GcmF1bWVuaSBzeW5kcm9tZSBzdWdnZXN0aW5nCiAgICAgICAgICAgICAgIyBnZXJtbGluZSBUUDUzIG11dGFudCBhcyBjYW5jZXIgcHJlZGlzcG9zaXRpb24gCiAgICAgICAgICAgICAgKHRwNTNfc2NvcmUgPiAwLjUgJiAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICBncmVwbCgiTGktRnJhdW1lbmkgc3luZHJvbWUiLGNhbmNlcl9wcmVkaXNwb3NpdGlvbnMpKSB8CiAgICAgICAgICAgICAgIyBjaGVjayBpZiB0cDUzIGluYWN0aXZhdGluZyBzY29yZSBmb3IgUk5BX1NlcSBpcyBncmVhdGVyIHRoYXQgMC41CiAgICAgICAgICAgICAgIyBhbmQgaGFzIDEgb3IgbW9yZSBTTlZ8IChDTlZ8U1YpIGxvc3MgaW4gVFA1MwogICAgICAgICAgICAgICh0cDUzX3Njb3JlID4gMC41ICYgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgKFNOVl9pbmRlbF9jb3VudHMgPj0gMSB8IChDTlZfbG9zc19jb3VudHMgPj0xIHwgU1ZfY291bnRzID49MSApKSkKICAgICAgICAgICkgfiAibG9zcyIsCiAgICAgICAgIyB3aGVuIGFjdGl2YXRpbmcgPT0gMQogICAgICAgIGFjdGl2YXRpbmcgPT0gMSB+ICJhY3RpdmF0ZWQiLAogICAgICAgICMgaWYgbm8gZXZpZGVuY2Ugc3VwcG9ydHMgVFA1MyBpbmF0aXZhdGlvbiBvciBhY3RpdmF0aW9uIAogICAgICAgIFRSVUUgfiAiT3RoZXIiCiAgICAgICkKICApCgpgYGAKCiMjIEV4cGxvcmUgZGlzdHJpYnV0aW9uIG9mIHRwNTNfYWx0ZXJlZCBzdGF0dXMgdnMgdHA1MyBpbmFjdGl2YXRpb24gc2NvcmVzCgpgYGB7ciwgb3V0LndpZHRoPSI1MCUifQpnZ3Bsb3QodHA1M19hbHRlcmF0aW9uc19zY29yZV93aWRlLCBhZXMoeCA9IGZhY3Rvcih0cDUzX2FsdGVyZWQpLCB5ID0gdHA1M19zY29yZSkpICsKICBnZW9tX3Zpb2xpbigpKwogIGdlb21faml0dGVyKGFscGhhID0gMC41LCB3aWR0aCA9IDAuMikgKwogIHN0YXRfY29tcGFyZV9tZWFucygpICsKICB0aGVtZV9idygpICsKICBnZ3RpdGxlKCJEaXN0cmlidXRpb24gb2Ygc2NvcmVzIGFjcm9zcyB0cDUzIGFsdGVyZWQgc3RhdHVzIikgKwogIHRoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gNjAsIGhqdXN0ID0gMSkpKwogIHhsYWIoInRwNTMgYWx0ZXJlZCBzdGF0dXMiKSAKCmBgYApUUDUzIG11dGFudHMgd2l0aCBnYWluLW9mLWZ1bmN0aW9uIChhY3RpdmF0aW5nKSBtdXRhbnRzIHByb21vdGUgdHVtb3JpZ2VuZXNpcyBhY2NvcmRpbmcgdG8gbGl0ZXJhdHVyZS4gW1JlZmVyZW5jZV0oaHR0cHM6Ly9wdWJtZWQubmNiaS5ubG0ubmloLmdvdi8xNzQxNzYyNy8pIGFuZCBbcmVmZXJlbmNlXShodHRwczovL3B1Ym1lZC5uY2JpLm5sbS5uaWguZ292LzI0Njc3NTc5LykgaGF2ZSBzaW1pbGFyIGRpc3RyaWJ1dGlvbiBvZiBjbGFzc2lmaWVyIHNjb3JlcyB0byB0aGF0IG9mIHBvdGVudGlhbCBsb3NzLW9mLWZ1bmN0aW9uIChtdWx0aS1hbGxlbGljKSBUcDUzIG11dGF0aW9ucy4gCgoiT3RoZXIiIGFubm90YXRlZCBzYW1wbGVzIGVpdGhlciBkb24ndCBoYXZlIGFueSBoaWdoLWNvbmZpZGVuY2UgbG9zcy9nYWluIFRQNTMgU05WcyBub3IgQ05WIGxvc3NlcyBPUiBETkEgc2FtcGxlIGlzIG5vdCBhdmFpbGFibGUuCgoKIyMjIEV4cHJlc3Npb24gcHJvZmlsZSBmb3IgYWN0aXZhdGluZyB2cyBsb3NzIFRQNTMgc3RhdHVzCmBgYHtyfQoKZXhwcl9jb21iaW5lZCA8LSByZWFkUkRTKGZpbGUucGF0aChkYXRhX2RpciwiZ2VuZS1jb3VudHMtcnNlbS1leHBlY3RlZF9jb3VudC1jb2xsYXBzZWQucmRzIikpCgpgYGAKCiMjIyMgKHN0cmFuZGVkKSBFeHByZXNzaW9uIHByb2ZpbGUgZm9yIGFjdGl2YXRpbmcgdnMgbG9zcyBUUDUzIHN0YXR1cwoKYGBge3J9CiMgc3Vic2V0IHRvIHN0cmFuZGVkIHNhbXBsZXMgCnN0cmFuZGVkX3NhbXBsZXMgPC0gY2xpbmljYWwgJT4lIAogIGZpbHRlcihzYW1wbGVfdHlwZSA9PSAiVHVtb3IiLAogICAgICAgICBleHBlcmltZW50YWxfc3RyYXRlZ3kgPT0gIlJOQS1TZXEiLAogICAgICAgICBjb2hvcnQgJWluJSBjKCJQQlRBIiwiR01LRiIsIlRBUkdFVCIpLAogICAgICAgICBSTkFfbGlicmFyeSA9PSAic3RyYW5kZWQiKSAlPiUKICBwdWxsKEtpZHNfRmlyc3RfQmlvc3BlY2ltZW5fSUQpICU+JSB1bmlxdWUoKQoKc3RyYW5kZWQgPC0gZXhwcl9jb21iaW5lZCAlPiUgc2VsZWN0KHN0cmFuZGVkX3NhbXBsZXMpCgojIHN1YnNldCB0byBUUDUzCnN1YnNldF9zdHJhbmRlZCA8LSB0KHN0cmFuZGVkKVssIlRQNTMiXQoKc3RyYW5kZWRfdHA1MyA8LSBhcy5kYXRhLmZyYW1lKHN1YnNldF9zdHJhbmRlZCkgJT4lCiAgZHBseXI6OnJlbmFtZShUUDUzPXN1YnNldF9zdHJhbmRlZCkgJT4lCiAgcm93bmFtZXNfdG9fY29sdW1uKCkgJT4lCiAgbGVmdF9qb2luKHRwNTNfYWx0ZXJhdGlvbnNfc2NvcmVfd2lkZSxieT1jKCJyb3duYW1lIj0iS2lkc19GaXJzdF9CaW9zcGVjaW1lbl9JRF9STkEiKSkgJT4lCiAgIyByZW1vdmUgcG9seWEgZnJvbSBtYXRjaGVkIHNjb3JlcwogIGZpbHRlcih0cDUzX2FsdGVyZWQgJWluJSBjKCJhY3RpdmF0ZWQiLCJsb3NzIikpCgpgYGAKClBsb3QgZGlzdHJpYnV0aW9uIG9mIFRQNTMgZ2VuZSBleHByZXNzaW9uIChzdHJhbmRlZCkKCmBgYHtyLCBvdXQud2lkdGg9IjUwJSJ9CgpnZ3Bsb3Qoc3RyYW5kZWRfdHA1MywgYWVzKHggPSBmYWN0b3IodHA1M19hbHRlcmVkKSwgeSA9IFRQNTMpKSArCiAgZ2VvbV92aW9saW4oKSsKICBnZW9tX2ppdHRlcihhbHBoYSA9IDAuNSwgd2lkdGggPSAwLjIpICsKICBzdGF0X2NvbXBhcmVfbWVhbnMoKSArCiAgdGhlbWVfYncoKSArCiAgZ2d0aXRsZSgiRGlzdHJpYnV0aW9uIG9mIFRQNTMgZXhwcmVzc2lvbiBhY3Jvc3MgdHA1MyBhbHRlcmVkIHN0YXR1cyIpICsKICB0aGVtZShheGlzLnRleHQueCA9IGVsZW1lbnRfdGV4dChhbmdsZSA9IDYwLCBoanVzdCA9IDEpKSsKICB4bGFiKCJ0cDUzIGFsdGVyZWQgc3RhdHVzIikgCgpgYGAKCiMjIyMgKHBvbHktYSkgRXhwcmVzc2lvbiBwcm9maWxlIGZvciBhY3RpdmF0aW5nIHZzIGxvc3MgVFA1MyBzdGF0dXMKCmBgYHtyfQoKIyBzdWJzZXQgdG8gcG9seWEgc2FtcGxlcyAKcG9seWFfc2FtcGxlcyA8LSBjbGluaWNhbCAlPiUgCiAgZmlsdGVyKHNhbXBsZV90eXBlID09ICJUdW1vciIsCiAgICAgICAgIGV4cGVyaW1lbnRhbF9zdHJhdGVneSA9PSAiUk5BLVNlcSIsCiAgICAgICAgIGNvaG9ydCAlaW4lIGMoIlBCVEEiLCJHTUtGIiwiVEFSR0VUIiksCiAgICAgICAgIFJOQV9saWJyYXJ5ID09ICJwb2x5LUEiKSAlPiUKICBwdWxsKEtpZHNfRmlyc3RfQmlvc3BlY2ltZW5fSUQpICU+JSB1bmlxdWUoKQoKcG9seWEgPC0gZXhwcl9jb21iaW5lZCAlPiUgc2VsZWN0KHBvbHlhX3NhbXBsZXMpCgojIHN1YnNldCB0byBUUDUzCnN1YnNldF9wb2x5YTwtIHQocG9seWEpWywiVFA1MyJdCgpwb2x5YV90cDUzPC0gYXMuZGF0YS5mcmFtZShzdWJzZXRfcG9seWEpICU+JQogIGRwbHlyOjpyZW5hbWUoVFA1Mz1zdWJzZXRfcG9seWEpICU+JQogIHJvd25hbWVzX3RvX2NvbHVtbigpICU+JQogIGxlZnRfam9pbih0cDUzX2FsdGVyYXRpb25zX3Njb3JlX3dpZGUsYnk9Yygicm93bmFtZSI9IktpZHNfRmlyc3RfQmlvc3BlY2ltZW5fSURfUk5BIikpICU+JQogIGZpbHRlcih0cDUzX2FsdGVyZWQgJWluJSBjKCJhY3RpdmF0ZWQiLCJsb3NzIikpCgpgYGAKClBsb3QgZGlzdHJpYnV0aW9uIG9mIFRQNTMgZ2VuZSBleHByZXNzaW9uIChwb2x5LWEpCgpgYGB7ciwgb3V0LndpZHRoPSI1MCUiIH0KCmdncGxvdChwb2x5YV90cDUzLCBhZXMoeCA9IGZhY3Rvcih0cDUzX2FsdGVyZWQpLCB5ID0gVFA1MykpICsKICBnZW9tX3Zpb2xpbigpKwogIGdlb21faml0dGVyKGFscGhhID0gMC41LCB3aWR0aCA9IDAuMikgKwogIHN0YXRfY29tcGFyZV9tZWFucygpICsKICB0aGVtZV9idygpICsKICBnZ3RpdGxlKCJEaXN0cmlidXRpb24gb2YgVFA1MyBleHByZXNzaW9uIGFjcm9zcyB0cDUzIGFsdGVyZWQgc3RhdHVzIikgKwogIHRoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gNjAsIGhqdXN0ID0gMSkpKwogIHhsYWIoInRwNTMgYWx0ZXJlZCBzdGF0dXMiKSAKCmBgYAoKCiMjIyMgKHBvbHktYSBzdHJhbmRlZCkgRXhwcmVzc2lvbiBwcm9maWxlIGZvciBhY3RpdmF0aW5nIHZzIGxvc3MgVFA1MyBzdGF0dXMKCmBgYHtyfQoKIyBzdWJzZXQgdG8gcG9seWEgc2FtcGxlcyAKcG9seWFfc3RyYW5kZWRfc2FtcGxlcyA8LSBjbGluaWNhbCAlPiUgCiAgZmlsdGVyKHNhbXBsZV90eXBlID09ICJUdW1vciIsCiAgICAgICAgIGV4cGVyaW1lbnRhbF9zdHJhdGVneSA9PSAiUk5BLVNlcSIsCiAgICAgICAgIGNvaG9ydCAlaW4lIGMoIlBCVEEiLCJHTUtGIiwiVEFSR0VUIiksCiAgICAgICAgIFJOQV9saWJyYXJ5ID09ICJwb2x5LUEgc3RyYW5kZWQiKSAlPiUKICBwdWxsKEtpZHNfRmlyc3RfQmlvc3BlY2ltZW5fSUQpICU+JSB1bmlxdWUoKQoKcG9seWFfc3RyYW5kZWQgPC0gZXhwcl9jb21iaW5lZCAlPiUgc2VsZWN0KHBvbHlhX3N0cmFuZGVkX3NhbXBsZXMpCgojIHN1YnNldCB0byBUUDUzCnN1YnNldF9wb2x5YV9zdHJhbmRlZDwtIHQocG9seWFfc3RyYW5kZWQpWywiVFA1MyJdCgpwb2x5YV9zdHJhbmRlZF90cDUzPC0gYXMuZGF0YS5mcmFtZShzdWJzZXRfcG9seWFfc3RyYW5kZWQpICU+JQogIGRwbHlyOjpyZW5hbWUoVFA1Mz1zdWJzZXRfcG9seWFfc3RyYW5kZWQpICU+JQogIHJvd25hbWVzX3RvX2NvbHVtbigpICU+JQogIGxlZnRfam9pbih0cDUzX2FsdGVyYXRpb25zX3Njb3JlX3dpZGUsYnk9Yygicm93bmFtZSI9IktpZHNfRmlyc3RfQmlvc3BlY2ltZW5fSURfUk5BIikpICU+JQogIGZpbHRlcih0cDUzX2FsdGVyZWQgJWluJSBjKCJhY3RpdmF0ZWQiLCJsb3NzIikpCgpgYGAKClBsb3QgZGlzdHJpYnV0aW9uIG9mIFRQNTMgZ2VuZSBleHByZXNzaW9uIChwb2x5LWEpCgpgYGB7ciwgb3V0LndpZHRoPSI1MCUiIH0KCmdncGxvdChwb2x5YV9zdHJhbmRlZF90cDUzLCBhZXMoeCA9IGZhY3Rvcih0cDUzX2FsdGVyZWQpLCB5ID0gVFA1MykpICsKICBnZW9tX3Zpb2xpbigpKwogIGdlb21faml0dGVyKGFscGhhID0gMC41LCB3aWR0aCA9IDAuMikgKwogIHN0YXRfY29tcGFyZV9tZWFucygpICsKICB0aGVtZV9idygpICsKICBnZ3RpdGxlKCJEaXN0cmlidXRpb24gb2YgVFA1MyBleHByZXNzaW9uIGFjcm9zcyB0cDUzIGFsdGVyZWQgc3RhdHVzIikgKwogIHRoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gNjAsIGhqdXN0ID0gMSkpKwogIHhsYWIoInRwNTMgYWx0ZXJlZCBzdGF0dXMiKSAKCmBgYAoKIyMjIENoZWNrIGlmIG90aGVyIGNhbmNlciBwcmVkaXNwb3NpdGlvbiBoYXZlIGhpZ2ggVFA1MyBpbmFjdGl2YXRpb24gc2NvcmVzICAKCmBgYHtyICxvdXQud2lkdGg9IjUwJSJ9CgpnZ3Bsb3QodHA1M19hbHRlcmF0aW9uc19zY29yZV93aWRlLCBhZXMoeCA9IGZhY3Rvcih0cDUzX2FsdGVyZWQpLCB5ID0gdHA1M19zY29yZSkpICsKICBnZW9tX3Zpb2xpbigpKwogIGdlb21faml0dGVyKGFscGhhID0gMC41LCB3aWR0aCA9IDAuMikgKwogIHRoZW1lX2J3KCkgKwogIGdndGl0bGUoIkRpc3RyaWJ1dGlvbiBvZiBzY29yZXMgYWNyb3NzIHRwNTMgYWx0ZXJlZCBzdGF0dXMiKSArCiAgdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA2MCwgaGp1c3QgPSAxKSkrCiAgeGxhYigidHA1MyBhbHRlcmVkIHN0YXR1cyIpICsKICBmYWNldF93cmFwKC5+Y2FuY2VyX3ByZWRpc3Bvc2l0aW9ucykKCmBgYApTb21lIE5GLTEgYW5kL29yIE90aGVyIGluaGVyaXRlZCBjb25kaXRpb25zIE5PUyBoYXZlIGhpZ2ggc2NvcmluZyBUUDUzIG11dGFudHMgYXMgd2VsbC4gCgpJbnRlcmVzdGluZyAyIGJzIGlkcyBhbm5vdGF0ZWQgd2l0aCBMaS1GcmF1bWVuaSBzeW5mcm9tZSBwcmUtZGlzcG9zaXRpb24gaGF2ZSB2ZXJ5IGxvdyB0cDUzIGNsYXNzaWZpZXIgc2NvcmVzLgoKYGBge3J9Cgp0cDUzX2FsdGVyYXRpb25zX3Njb3JlX3dpZGUgJT4lCiAgZmlsdGVyKGNhbmNlcl9wcmVkaXNwb3NpdGlvbnMgPT0gIkxpLUZyYXVtZW5pIHN5bmRyb21lIiwKICAgICAgICAgdHA1M19zY29yZSA8IDAuNSkgCiAgCmBgYAoKCiMjIFNhdmUgZmlsZSAKCkFkZGluZyBETkEgYW5kIFJOQSBiaW9zcGVjaW1lbiBpZHMgdG8gYHRwNTNfYWx0ZXJhdGlvbnNfc2NvcmVfd2lkZWAgYW5kIHNhdmluZwoKYGBge3J9Cgp0cDUzX2FsdGVyYXRpb25zX3Njb3JlX3dpZGUgJT4lIAogIHdyaXRlX3RzdihmaWxlLnBhdGgocmVzdWx0c19kaXIsInRwNTNfYWx0ZXJlZF9zdGF0dXMudHN2IikpCgpgYGAK